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A PROCEDURE FOR OBTAINING STIFFNESSES 
AND MASSES OF A STRUCTURE FROM VIBRATION 
MODES AND SUBSTRUCTURE STATIC TEST DATA 

Harold Edighoffer* 

Langley Research Center 

SUMMARY 

The purpose of this investigation was to develop a component mode 
desynthesis procedure for determining the unknown vibration character- 
istics of a structural component (i. e. , a launch vehicle) given the vibra- 
tion characteristics of a structural system composed of that component 
combined with a known one (i. e. , a payload). The desynthesis procedure 
developed does not require the vibration characteristics of the payload 
component but requires that at least one component static test be per- 
formed. This data is used in conjunction with the system measured 
frequencies and mode shapes to obtain the vibration characteristics of 
each component. The static test data is a catalyst in the de synthesis 
procedure and could be performed on any component or on the entire 
system. This method has direct application to determining the flight 
dynamics of an empty launch vehicle from measurements made on a 
vehicle/payload combination in conjunction with a static test on the payload. 
This procedure is applicable to a multicomponent system. 
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L INTRODUCTION 


This model desynthesis procedure consists of the following steps. 
First, the coupled degrees of freedom for the structural system 
is established. Secondly, a component static test is performed. The 
unknown system stiffness and mass matrices are determined by parts, 
a row at a time, using the system measured eigenvalues and eigenvectors 
with the static test data. The stiffness and mass matrices of each com- 
ponent and the interface springs are determined from the system matrices. 
Finally, the component mass and stiffness matrices of the components 
are used in an eigenvalue analysis to determine the component dynamic 
responses. The method is demonstrated by the evaluation of a two 
beam, twelve degree-of-freedom system. 
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SYMBOLS 


K or K. . 

P iJP 

K or K-. 
ij 


M or M 






ut or Uit 


u or u 


ir 


2 

u) 

r 


0 ) 





Payload stiffness 

System stiffness 

System mass 

Static test force 

Static test displacement 

Eigenvector, i = degree of freedom, 
r = mode number 

Eigenvalue 


6^or 6^ 


Measured static test 


“t' ®i Pi 
F or P 

T 

b. 


Measured static test 


Measured static test 


Measured static test 


Column matrix 


superscript or subscript 
r Mode number 


displacement 

rotations 

force 

torque 
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III. DESYNTHESIS PROCEDURE 


3. I Determine the Coupled Degrees of Freedom 

The first step in this desynthesis procedure is to determine all 
of the structural coupling terms in the mass and stiffness matrices. If 
the structure has stiff bulkheads so that they respond as rigid bodies, 
the first and third bulkhead are uncoupled in the mass and stiffness 
matrix. Furthermore, if the motion of one bulkhead in a given direction 
results in a net force of zero on an adjacent fixed bulkhead these two 
degrees of freedom are uncoupled. The unknown stiffness matrix with 
an X representing the coupled degrees of freedom is shown in equation (1) 
for a 12 degree of freedom system. 

Coupled terms in a twelve degree-of-freedom system . - The two 
beam model is shown below. The model properties are shown in Figure 1. 
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The coupled degrees of freedom in K matrix are; 



For this structure there is no coupling between degree of freedom 1 and 
5, 6, 8, 9, 10, 11 and 12 because when the degree of freedom 1 is 
given a unit displacement with all others restrained, there is no reaction 
force at these uncoupled degrees of freedom. By walking through each 
degree of freedom in a similar manner for degree of freedom 2, 3, etc. , 
the complete system unknown coupled K term locations are determined. 
There is no interface spring coupling between 1-8, 2-7, 5-12 and 6-11 
because of knowledge of the interface spring design. There is no coupling 
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between 3-4 and 9-10 because of the mid point location with uniform stiff- 
ness along the length. The mass matrix will have the identical coupled 
locations except the 1-7, 2-8, 5-11 and 6-12 couplings which are omitted 
because there is no mass coupling across the interface. 

3. 2 Component Static Test 

The second step is to perform a static test on one of the com- 
ponents. The payload component will normally have fewer degrees of 
freedom then the launch vehicle component and would probably be more 
economical to test. Alternatively, the entire system could be statically 
tested to satisfy this requirement. 

The component static test data can be expressed in matrix form as: 

[Kp] (2) 

where: 

K = unknown payload stiffness matrix 

Uj. = static test displacements 

F j. = static test forces (applied and reaction) 


For the twelve degree-of-freedom system equation (2) becomes: 
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and zero but were included in this analysis to demonstrate 

that uncoupled degree-of-freedom combinations can be overlooked to 
begin with and subsequently identified by the method. 

The terms with the subscript p are the payload stiffness terms that 
are different than the system terms in equation (1) because of the coupling 
springs. The system terms in equation (1) are related to the payload 
terms of equation (3) and the interface stiffness terms by 


1 = K 

lip 

- Ki7 

^22 = ^22p 

- ^28 

1^55 = ^55p 

- ^5. 11 

^66 = K66p 

- ^6. 12 


Static test for payload of twelve degree-of-freedom system . - The 
payload component is simply supported at each end and loaded in bending 
as shown below. Other tests such as a cantilever beam test would be just 
as acceptable. 



6 


With a test force of F, the payload component has a measured 6j- 
deflection at the mid-point and a measured at rotation at each end. The 
test displacements and forces of equation (3) becomes: 



Equation P) is revised by including the equalities of equation (4) to 
express all stiffness terms in system notation. 
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Imposing the condition that K^j = Kjx fori Equation (6) 

can be rearranged row by row, so that the unknown stiffnesses are 


contained in vectors (column matrices). 
1st row 



(7a) 


(7b) 


(7c) 
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4th row 



this point it is necessary to decide which row of the system K matrix is 

to be solved first. Because of the zero static test forces in equation (5), 

the 2nd, 4th and 6th rows equations (7b), (7d), and(7f) are eliminated, since the 


9 



purpose of the static test is to provide finite values on the right hand side of 
a set of simultaneous equations so as to obtain stiffnesses (and also masses). 
For this example, the first row (equation 7a) was selected for the first set of 
simultaneous equations. 

3.3 System Measured Eigenvalues and Eigenvectors 

The system stiffness and mass matrices are related to each measured 
eigenvalue and its associated eigenvector of the vehicle/payload system as 



where 

K = unknown vehicle/ payload stiffness matrix 

M = unknown vehicle /payload mass matrix 
2 

= the rth known system eigenvalue 

(r) 

u = the rth set of eigenvectors 

Revise equation (8) for the twelve degree-of-freedom system. - For 
the example two-beam system, equation 8 for the rth eigenvalue becomes 
equation (9) on the next page. 
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"2r 

^3r 

“4r 

^5r 

'^6r 

U9r 

“lOr 

'^llr 

V_“l2rJ 


ru_■^ 


u 


Ir 


2r 

'iSr 

U4r 


> -u, 


(r) 




U, 


5r _ 


u 


6r 


“7r 

“8r 

U9r 

'^lOr 

“Hr 

Cl2r^ 


>= ^ 


^0^ 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 






(9) 


where the system K matrix of unknown coupled terms are identical to 


equation (1) with the 3-4 and 9-10 terms added. 


^11 

«12 

^13 

^14 

0 

0 

^17 

0 

0 

0 

0 

0 

Kzi 

«22 

K23 

^24 

0 

0 

0 

Kz8 

0 

0 

0 

0 

K31 

'^32 

K33 

^34 

K35 

^36 

0 

0 

0 

0 

0 

0 

^41 

“^42 

K43 

K44 

K45 

^46 

0 

0 

0 

0 

0 

0 

0 

0 

^53 

^54 

K55 

^56 

0 

0 

0 

0 

K5.11 

0 

0 

0 

^63 

^64 

^65 

K66 

0 

0 

0 

0 

0 

^,12 

^71 

0 

0 

0 

0 

0 

K77 

^8 

K.,„ 

79 

Vo 

0 

0 

0 

K82 

0 

0 

0 

0 

^87 

^88 

^89 

0 

f 

00 

0 

0 

0 

0 

0 

0 

0 

0 

K97 

K98 

K99 

^,10 

K9.11 

K9.12 

0 

0 

0 

0 

0 

0 

^10,7 

^10,8 

^10,9 

^laio 

^10.11 

^10.12 

0 

0 

0 

0 

Kii ,5 

0 

0 

0 

Kii,9 

^11,10 

K11.11 

Kii,12 

0 

0 

0 

0 

0 

^12,6 

0 

0 

^12,9 

^12,10 

^12.11 

^12.12 


11 



The system mass matrix is identical except the 1-7, 2-8, 5-11 and 6-12 
terms are zero. 

Imposing the symmetry condition that s Kj^ and M^j s for i j. 

row of equation (9) is rearranged to be 
1 st row 


^Ir “2r "3r “4r "7r ^Ir ^Zr ^3r ^4r j 


where 


W. = -U) u- 
”ir r ir 


K 


K 


K 


K 


11 


12 


13 


14 


K 


17 


M 


11 


M 


12 


M 


13 


Mi4 


= 0 


(10a) 


2nd row 


l"lr “zr “3r “4r “8r ""ir '^2r ^3r j 


where 


Wir = -<“r «^ir 


^12 


K 


K 


22 


23 




K 


24 


K 


28 


M 


12 


M 


22 


M 


23 


^24 


= 0 


(10b) 


Each 
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3rd row 



(10c) 
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5th row 


I ^3r «4r "Ur ^ 3 ,. 


where 


Wir = -'"r "ir 


^ K35 >1 


K 


45 


K 


55 




K 


56 


K 


> =0 (lOe) 


5,11 


M 


35 


M 


45 


M 


55 


M 


56 


6th row 


|"3r '^4r "5r "6r '^izr ^3r ^4r ^5r ^6r | 


where 


Wir = -«Jr "ir 


K 


36 


K 46 
K 56 

■^66 

<^ 6,12 
M 36 
M 


> =0 (lOf) 


"46 


M 


56 


M 


66 
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7th row 



8th row 



(lOh) 


15 



9th row 



16 



1 1th row 
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3.4 


Formulate Simultaneous Equations for Each Row of K and M Matrix 
Each of the equations (10a) thru (lOj) can be expanded by letting 


r = 1, 2, . . . , m, to obtain a set of simultaneous equations that could 
be solved from the measured system eigenvalues and eigenvectors except 
that all the equations have a zero on the right hand side. To circumvent 
this problem one static test data equation (equation 7) is combined with 
the number of dynamic response equations (equation 10) required to equal 
the number of unknownK and M terms. Obviously, the right hand side of 
the static test equation (7) used, must not equal zero. For this example 
12 degree-of-freedom system, equation (7a) for the first row was combined 
with (lOa). Equation (lOa) has 9 unknown system K and M terms. By com- 
bining the first 8 eigenvalues (r = 1, 2, 8 in equation 10a) 

with the first row test equation (equation 7a), the 10 values are solved 
s imultaneou sly. 

Combine Equation (7a) and (10a) for the 1st row simultaneous 
equations . - This is assembled in matrix form to solve for the first row 
of the mass and stiffness matrix as shown in equation (11). 
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1st Row Simultaneous Equations 



^21 

“31 

“41 

“71 

0 

0 

0 

0 


r >v 
^11 


0 

"12 

“22 

^32 

“42 

“72 

0 

0 

0 

0 


Ki2 


0 

^13 

“23 

“33 

“43 

“73 

^13 

^23 

^33 

^43 


Ki3 


0 

ui 4 

“24 

“34 

“44 

“74 

Wi4 

W24 

W34 

W44 

i 

^14 


0 

“15 

“25 

“35 

“45 

“75 

^15 

^25 

W35 

^45 


Ki7 


0 

^16 

“26 

^^36 

“46 

“76 

^16 

^26 

^36 

^46 


Mil 


0 

“17 

“27 

“37 

“47 

“77 

^^17 

W27 

^37 

W47 


^12 


0 

“18 

“28 

“38 

“48 

“78 

^18 

^28 

W38 

^48 


^13 


0 

^It 

“2t 

“3t 

“4t 

“It 

0 

0 

0 

0 




Fit 

'w ) 


where 


W = 
ir 


CD 


ir 


The first two rows are rigid body modes and the last row is from the 
static test data, so 6 flexible system test eigenvalues and associated eigen- 
vectors are required. An alternate approach would be to use 8 test modes 
with a static test equation and eliminate the rigid body modes. 

The rigid body modes are constructed. - The first rigid body mode 
is pure translation and the 2nd is rotation about the mid points of both 
components. These rigid body eigenvalues are zero and the eigenvectors 
are obtained directly from the geometry of the structure. Other rigid body 
displacements could be used. 
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Combine equation (7b) and (10b) for the second row simultaneous 
equations , - In a manner similar to the derivation of equation (11), the 
K 2 « and terms are solved by transposing the known Kj 2 ^IZ 
that were solved in equation (11). There are 7 unknowns remaining. 


Znd Row Simultaneous Equations 


“23 

'^33 

"43 

"83 

^23 

W33 

W4T 


^22 




"13 

W 

13 

"24 

"34 

"44 

"84 

^24 

W34 

W44 


^23 


0 


"14 

^14 

"25 

“35 

"45 

"85 

^25 

^35 

^45 

i 

^24 

i 

0 

>- 

"15 

^5 

"26 

"36 

"46 

"86 

^26 

W 

'^36 

^46 


^28 


0 


"I6 

^16 

*^27 

"37 

"47 

"87 

^27 

W 

37 

W47 


^22 


0 


"17 

^17 

"28 

"38 

"48 

«88 

^28 

W38 

W48 


M2 3 


0 


"18 

^18 

"2t 

"3t 

"4t 

"2t 

0 

0 

0 


^24 


®2t 


"it 

0 


where 



(12) 


W. = -u) 
ir r 


u 


The two rigid body modes were not used. Other options are to use one 
rigid body mode and eliminate the static test equation or use 7 flexible 
modes. 


Combine equation (7c) and (10c) for the third row simultaneous 
equations . - The ^13 terms were solved in equation (11) and the 

K 23 and ^23 terms were solved in equation ( 12 ) and are transposed to the 
right hand side. There are 8 unknowns remaining. 
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3rd Row Simultaneous Equations 


"32 "42 


U 33 U 43 


"34 "44 


U.p 

35 45 


""36 "46 


"37 "47 


"38 "48 


_"3t "4t 

where 


"52 "62 " u u u 

"53 "63 "^33 "^43 "^53 "^63 

"54 “64 ^^34 '^44 ’^54 ^^64| 

"55 "65 '^35 '^45 '^ S5 '^65 

“56 “66 ^^36 '^46 ^^56 ^^66 

"57 “67 '^37 ^^47 '^57 ^7 

“58 “68 "^38 "^48 '*^58 ^68 




^33 


K 




34 
*^35 

"36 ) 
M 


33 


M 


34 


M 


35 


"5t " 

6t 0 


0 0 

»J 

1 

fo' 


r 

"12 

“22 

0 


: ° 


"13 

"23 

^13 

w 

^23 

1 0 


"l4 

“24 

^14 

W 

24 

' 0 

^ - 

"15 

"25 

^15 

Wz5 

0 


"16 

""26 

^16 

^26 

0 


"17 

"27 

W 

17 

^27 

0 


"18 

"28 

W18 

^28 

..^3t> 


"it 

^2t 

0 

0 


M 


36 




^13 


K 


23 


M 


13 


V. 23y 


('’3) 


(13) 


W. = - oj u, 
ir r ir 


One rigid body mode and the static test was used. Other options would be 
to use the two rigid body modes and eliminate the static test equation or 
use 8 flexible modes. 
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Combine equation (7d) and (lOd) for the fourth row simultaneous 


equations ^ - and from equation (13) are transposed to the right 

hand side. There are 6 unknowns remaining. 


4th Row Simultaneous Equations 


^43 

^53 

^"63 

^43 

^53 

W " 
63 


44 

^44 

“54 

“64 

W44 

^54 

^64 

i 

K45 

'45 

“55 

“65 

W45 

W55 

^65 


K46 

46 

“56 

“66 

^46 

^56 

^66 


M44 

■47 

“57 

“67 

^47 

^57 

^67 


M45 

•48 

“58 

“68 

^48 

W58 

^68 





where 



“13 

“23 

"33 

^3 

^23 

^33"' 


Ki 4 

“14 

“24 

"34 

^14 

^24 

^34 


K24 

ui5 

"25 

"35 

Wi5 

^25 

W35 


^34 

^^16 

“26 

"36 

^16 

^26 

W 

36 

A 

Mj4 

“17 

“27 

"37 

Wi7 

^27 

W37 


^24 


“28 

“38 

^18 

^28 

^38 


< 34J 


W. 


ir 




u. 

ir 


(14) 


The first 6 flexible modes were used. Other options would be to include 
the static test equation and/or the rigid body modes. 
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Combine equation (7e) andflOe) for the fifth row simultaneou s 


equations ^ - ^ 25 * ^ 25 ' ^35' ^35' ^45' ^45 pre- 

viously been solved in equation ( 11 ) thru ( 14 ) and are transposed to the right 
hand side. There are 5 unknowns remaining. 


5 th Row Simultaneous Equations 


*^53 

"63 

"ll,3 

"^53 

^63“ 


"K55' 


"33 

“43 

33 

^3' 



*^54 

“64 

"ll,4 

"'54 

^64 

i 

^56 

►= - 

“34 

“44 

W 

‘ 34 

^44 

< 

«45 


“65 


w 

55 

W 

65 


^,11 


"35 

“45 

W 

35 

^^45 


M35 

"56 

“66 

“ll,6 

W 

56 

W 

66 


M55 


“36 

"46 

W 

36 

w 

46 


M45 


u 

u 


W 




u 


W 

W 



57 

67 

11,7 

57 

67^ 




- 37 

47 

37 

47 




where 


W. = -u> 
ir r 


u 


ir 


The first 5 flexible modes were used. Singularity problems were en- 
countered when both rigid body modes and the test equation was uee,d 
with 2 flexible modes. 
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Combined equation (7f) and (lOf) for the sixth row simultaneous 
equations . - ^56' ^36' ^46 ^56 previously been 

solved in equations (13) thru (15) and are transposed to the right hand side. 
There are 3 unknowns remaining. 


6th Row Simultaneous Equations 


"63 

'"12,3 

^63' 


fv . > 
^66 


"33 

"43 

"53 

W33 

W43 

^53" 

"64 

'"l2,4 

^64 

< 

^6.12 

> = - 

“34 

"44 

"54 

^34 

W44 

W54 

"65 

'"l2, 5 

^65 


^^66^ 


“35 

*"45 

"55 

^35 

W45 

^55_ 


where 




K 


36 


K 


Ks6) 

M, 


36 


^46 

V, 


w.^ = 


CU 


U 


ir 


The first 3 flexible modes were used. One rigid body mode or the static 
test equation could be used with two of the flexible modes. The number 
of flexible modes used should be greater than the sum of the rigid body 
modes to eliminate singularity problems. 
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Expand equation (lOg) for the seventh row simultaneous equation s. - 
has previously been solved in equation (1 1) and is transposed to the right 


hand side. There are 8 unknowns remaining. 
7th Row Simultaneous Equations 


“71 

“81 

“91 

“10, 1 

0 

0 

0 
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"*^ 77 ^ 


“11 

“72 

“82 

“92 

“10,2 
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W73 
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W93 
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U74 
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74 

^84 
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^10.4 
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u 

14 

“75 

“85 

“95 

“10,5 

W75 

^85 

W95 

^10, 5 


M77 


“l5 

“76 

“86 

“96 

“10,6 

^76 
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Wl0,6 


M78 


“16 

“77 

“87 

“97 

“10,7 

W77 

^87 

W 

97 

^10,7 


M79 


“17 

00 

“88 

“98 

“10, 8 

W78 

^88 

W 

98 

^10,8 




^“18^ 


K 


17 


where 


W. = - u. 
ir r ir 


(17) 


This IS the first row of K and M terms in the launch vehicle component that 
was not tested statically. Six measured system eigenvalues and eigenvectors 
were used with the 2 rigid body modes. Eight flexible modes could be used 
with no rigid body modes. 
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Expand equation (1 Oh) for the eighth row simultaneous equations . 
^28' K^g and have previously been solved in equations (12) and 
(17) and are transposed to the right hand side. There are 6 unknowns 
remaining. 


8th Row Simultaneous Equations 
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“93 

“lO, 3 

^83 

W 93 

^ 10 , 3 " 
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^10,5 
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'"76 
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“27 

“77 
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^98 

“l0,7 
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W 98 

Wio,8 


^8,10^ 


“28 

“78 
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where 


W. 


ir 


2 

r ir 


The first 6 flexible modes were used. The 2 rigid body modes could have 
been used. 
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Expand equation (lOi)for the ninth row simultaneous equations. - 
Kyg, Kgg, ^79 ^89 previously been solved in equations (17) and 

(18)and are transposed to the right hand side. There are 8 unknowns 
remaining. 


9 th Row Simultaneous Equations 
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“11,1 

“12, 1 

0 

0 

0 

0 


K99 

^92 

‘^10,2 

*^11,2 

^^12, 2 

0 

0 

0 

0 


^9,10 

^93 

^10, 3 

^11,3 

^12,3 

^9,3 

^10,3 

^11,3 

^12, 3 


^9,11 

U94 

’^10,4 

‘^11,4 

^^12, 4 

^9,4 

^10,4 

^11,4 

^12,4 


S, 12 

'"95 

'^lO, 5 

'^11,5 

'"12, 5 

W9.5 

^10,5 

^11, 5 

^12,5 


M99 

^^96 

‘^10,6 

'"11, 6 

*^12,6 

^9,6 

^10,6 

^11,6 

^12,6 


^9,10 

U97 

1 

‘^10, 7 

'"11,7 

^^12, 7 

^9,7 

Wio,7 

W11.7 

^12, 7 


^9,11 

^^98 ‘^10,8 
where 

'"11,8 

'"12, 8 

^9,8 

^10. 8 

Wii,8 

^12,8 

Ml 

• 


^^9,12^ 






(19) 




U 


u 


u 


ir 


i) ^ 

r ir 


71 

"81 

0 

0" 


( Vi 

K79 

72 

"82 

0 

0 


K89 

73 

""83 

^73 

^83 

M79 

74 

'"84 

^74 

00 


00 

^ 

75 

'"85 

W 

75 

^85 



76 

"86 

W76 

^86 



77 

"87 

W77 

00 



78 

"88 

W78 

00 ^ 
00 



u. 







The 2 rigid body modes and the first 6 flexible modes were used. Another 
option would be to use 8 flexible modes. 
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Expand eq uation gOj) for the tenth row simultaneous enuations. - 
K?, 10» jQ, 10 , iQ, M 3 ^ 10 and jo have been solved in 
equations (17), (18) and (19) and are transposed to the right hand side. There 
are 6 unknowns remaining. 

10th Row Simultaneous Equations 


"10,3 

"11,3 

"12,3 

^10,3 

^11,3 

^12,3" 


(Sc ^ 

^10, 10 

"10,4 

“ll,4 

"12,4 

Wi0,4 

^11,4 

^12,4 


^10, 11 

“l0,5 

"ll,5 

"l2,5 

^10,5 

^11,5 

^12,5 

< 

^10, 12 

"10,6 

"ll,6 

"12,6 

^10,6 

^11,6 

^12,6 


^10, 10 

"lO, 7 

^^11,7 

“l2, 7 

^10,7 

^11,7 

^12,7 


^10, 11 

"lO, 8 

^^11,8 

"l2, 8 

W 

10,8 

^11,8 

^12, 8^ 


J^10,12^ 

where 

f"73 

"83 

"93 ^73 

; ^83 

W93I 

fK; 

lol 


= {»io} 


(20) 


{bio)= - 


U74 Ug4 Ug4 W74 Wg4 Wg4 

"75 "85 '^95 ^85 ^95 

"76 «86 "96 ^76 ^86 ^95 
"77 "87 "97 W77 Wg^ Wg7 


"78 "88 "98 ^78 ^88 


W 


98 


K 


8, 10 


K 


9, 10 


M 


M 


7, 10 

8, 10 


M 

V. 


9, 10 


W; = - 0) U; 
ir r ir 


The first 6 flexible modes were used. Another option would be to use 


2 rigid body modes with 4 flexible modes. 



Expand equation (10k) for the eleventh row simultaneous equations^ 
K 5 , li; Kg jj; Kjq jj; Mg jj and have been solved in equations 

(15), (19) and (20) and are transposed to the right hand side. There are 4 
unknowns remaining. 

11th Row Simultaneous Equations 


"^11, 3 

^12,3 

'^11,3 

^12,3 


''" 11 . ill 


‘^53 

“93 

“10,3 

Wg3 

^10.3 


^TC 

^5.11 

"11,4 

^12,4 

^11.4 

^12,4 


^11,12 

H- 

^^54 

'^94 

“10,4 


W10.4 


Hll 


^12.5 

^11,5 

^12,5 


^11, 11 


“55 


'^10,5 

"^95 

W 

10,5 

< 

K 

10,11 

""11,6 

‘"12,6 

^11,6 

^12, 6_ 




_‘^56 


“10,6 

% 

W 

10,^ 


^9,11 

who re 










^^10,11, 


( 21 ) 


W = - u) u. 
ir r ir 


The first 4 flexible modes we^e used. 
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Expand equation MOml for the twelfth row simultaneous equations. 


^ 6 , 12 » ^ 9 , 12 ’ ^ 10 , 12 ' ^ 11 , 12 ' ^ 9 , 12 ’ ^ 10, 12 12 

solved in equations (16), (19), (20) and (21) and are transposed to the right hand 
side. There are two unknowns remaining. 


12th Row Simultaneous Equations 


"12,2 Wj 2^2 


“62 “92 “10,2 “ll,2 '^10,2 

"12,6 ^12,6 

1.^12, 12 j 

“66 “% “10.6 “11,6 '^lO.O '^11,6. 


where 


K 


6,12 


K 


9,12 


K 


K 


M 


M 


10,12 

11,12 

9,12 

10,12 


(22) 


k.^ll,12j 


W. = -U) u 
ir r ir 


The second rigid body mode and the fourth flexible moae were used. 
The fourth flexible mode was used because of a high rotational degree 
of freedom IZ dynamic response. 


30 



IV. ANALYSIS AND DISCUSSION 


The twelve degree-of-f reedom two beam model shown in Section 3. 1 
and Figure 1 was evaluated using the simultaneous equations (11) thru (2Z)* 

The simple beam static test shown in Section 3.2 was used in this evaluation. 

The desynthesis computer code SIMUL3 was written. - The program 
uses the NASA Langley library program ITIMP to solve the simultaneous 
equations. The program input is the measured eigenvalues and eigenvectors, 
the rigid body eigenvectors, and the static test displacements and forces. 

The program stashes the data in matrix form into equation (11) and solves 
for the 9 values in the first row of the K and M matrices. The program then 
uses this output to evaluate the second row using equation (12) and continues 
through equation ^2) when all of the system stiffness and mass properties 
are determined. 

System dynamic data was simulated by a NASTRAN analysis . - To 
represent the first 6 measured system flexible eigenvalues and eigenvectors, 
the NASTRAN program was used with the eigenvalues normalized to unity 
(MAX and the NASTRAN EIGR card). The NASTRAN data for the 6 flexible modes 
and 2 rigid body modes are tabulated in Table 1. 

The payload static test data was simulated by a NASTRAN anaylsis. - A 
simple beam unit load (F= 1. 0) was applied at the mid point and the simulated test 
values of the end slope ^rt) center deflection (6t) of equation (5) were 

calculated using the properties in Figure 1. ' 

Four different evaluations using the computer code SIMUL3 were made . - 

Analys IS No, I . - The use of static test data with no measurement errors 
was evaluated with the system dynamic properties from Table 1. The stiff- 
ness and mass terms were compured with the correct values (from NASTRAN). 

The comparison is shown in Table 2 as a percentage change in mass or 
stiffness. No problems were encountered until the evaluation of 
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(equation 16) and Mj2, 

12 (equation 22). The value was only 84. 1% 

of the NASTRAN value but the original ^ 12 , 12 

was 1366% higher than 

the NASTRAN value so 

different options were used to evaluate Kj 2 12 

^12, 12 ^sing equation (22). 



% Change 

f D®.^ni_he^l^ ^ jqqI 
|_ NASTRAN J 

Inputs 

Ki2, 12 

^12. 12 

1 . 1st flex, mode 
2nd flex, mode 

99.6 

1366.2 

2. 2nd flex, mode 
3rd flex, mode 

99. 0 

-117.4 a 

3. 2nd rigid body 
1st flex, mode 

100.4 

-HIGH ^ 

4. 2nd rigid body 
2nd flex, mode 

100.4 

184. 1 

5. 2nd rigid body 
4th flex, mode 

100.4 

103.6 


The final values used were from the fifth input combination of the 2nd rigid body 
mode and the 4th flexible mode. The variation depends on what two modes 
are used and could be due to many variables. The use of the 4th flexible 
mode could have improved the results because this mode has a higher rota- 
tional dynamic response then the 1st, 2nd and 3rd flexible modes. Another 
possibility is that numerical problems because of small difference of large 
numbers in equation (22). The most likely cause could be the order in which the 
simultaneous equations were solved (equation 11 thru 22). The and 

Mil 11 terms have very large values relative to the Mj^2, 12 rotational 
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a. Negative values indicate the Desynthesis analysis value had the wrong sign. 



term indicating that the ^12 12 solved last.^ If the order of 

solving the rows in the K and M matrix are changed so that the Mgg 

or 11 last equation solved, this Kj2^ 12 ^md 12 discrepancy 

will probably be eliminated. Myy is 39 times larger than Mj2^ 12 
IS 78 times larger than Mj2, 12* Examination of Table 2 indicates that the 
same problem exists for the terms. The rearrangement of the order 

of analysis to solve for the before M 55 will probably eliminate the 

discrepancy. 

Even though this problem had not been eliminated, the K and M 
values of the payload component were evaluated for their dynamic response 
using SCAMP and compared with the dynamic response using the K and M 
values obtained from NASTRAN. The expected poor comparison is shown 
in Table 6 . 

Analysis No. 2. - The static test slope (Q^ |.) measurement at each 
end was increased 1 % and the displacement at the center was maintained at 
the correct value. The resulting K and M terms are compared with the 
results of Analysis No. 1 in Table 3. All of the K and M terms in the first 
row from equation (I Dincreased to 103. 1 % of the Analysis No. 1 results. 

The values in the second row varied from 84.2 to 106. 7%. The values in 
the third row varied from 32. 2% to 123. 3%. The error increased as the 
analysis continued through the payload component. However, when row 7, 
the first row in the vehicle component, was evaluated using equation (17), 
the error decreased back to 103. 1%, the value identical to row 1. In general, 
all of the stiffness terms become larger and the mass terms become smaller. 
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b. See Appendix for verification by subsequent analysis. 



The erratic changes in ^IZ, IZ ^IZ, 12 probably caused 

by the same problems that was encountered and discussed in Analysis No. 1. 
The and actually changed sign. The K or M matrix of the vehicle 

component was entered into SCAMP to determine the dynamic response 
and it was determined that the mass matrix was not positive definite and 
no eigenvalue solution was obtained. 

Analysis No. 3. - The static test displacement (6 measurement 
at the mid point was increased 1% and the slopes at the end points were 
maintained at the correct value. The resulting K and M terms are compared 
with the results of Analysis No. 1 in Table 4. All of the K and M terms in 
row 1 decreased to 96. Z% of the Analysis No. 1 values. After the first row 
the K terms tended to decrease and the M terms tended to increase in value. 
In the vehicle component the 7th row values were again the same as the row 1 
values at 96. Z%. No eigenvalue solution was attempted. 

Analysis No. 4. - The combined static test measurement error of 
Analysis No. Z and Analysis No. 3 was used. This included an increase in 
the slope at each end of 1% and an increase in displacement at the center 1%. 
The resulting K and M terms are compared with the results of Analysis No. 1 
in Table 5. As can be seen every value was 99. 0% of that in Analysis No. 1. 
Putting these vehicle component M and K terms into SCAMP resulted in 
identical eigenvalues to that of Analysis No. 1 as shown in Table 6. If the 
static test measurements are changed so that all readings are high or low 
the same percentage, the desynthesis analysis will give the correct vehicle 
component dynamic response. 
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V. EVALUATION OF LANGLEY DYNAMIC RESEARCH MODEL 


A preliminary evaluation was made to investigate the use of the Langley 
Dynamic Research Model in this desynthesis procedure. There was some 
system free-free dynamic response test data and some static test data avail- 
able on the payload. This model has a vehicle carrier and a payload com- 
ponent. The payload component shown in Figure 2 has 48 degrees of freedom. 

The payload and carrier together have 288 degrees of freedom. There are 
8 triangular bulkheads that are fairly rigid, connected by three 
longerons at the three corners. There are additional diagonal 
elements not shown to give the structure the shear and torsional capability. 

Each bulkhead can be considered to move as a rigid body with six degrees of 
freedom each or a total of 48. The degrees of freedom of the interconnecting 
springs between the payload and the carrier are shown in Figure 2. To 
gain some insight into the type of static test required, the uncoupled stiff- 
ness (K) terms are shown for the payload as an X at each coupled location 
in Figure 3. The first 6 degrees of freedom represent the first bulkhead in 
the six system directions (3 translations and 3 rotations). The second set of 
6 degrees of freedom (7 thru 12) represent bulkhead 2, etc. By studying 
Figure 3, row by row, it can be seen that the 1st bulkhead is coupled to the 
2nd bulkhead (rows 1 thru 6). The second bulkhead is coupled to the 1st and 
3rd bulkhead (rows 7 thru 12). Each bulkhead is coupled to the bulkheads 
adjacent to them. Also, some of the 6 degrees of freedom within each bulk- 
head are uncoupled. The diagonals between bulkheads (not shown) in figure 2 are 
arranged so that if the first bulkhead has a rigid bydy motion in the 3 direction with 
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the second bulkhead fixed, the net load in the 1 direction on bulkhead 1 is 
zero. Therefore, there are no 1-3 coupling terms in Figure 3. When the 
first bulkhead is moved in the 2 direction there is a net force in the 1 direc- 
tion. Therefore, the 1-2 coupling terms are present. Rows 3, 25, 26, 27, 
28, 29, and 30 have interface spring coupling terms not shown in Figure 3. 

The static test equation . - From equation (2) and equation (3), the 
first SIX rows of a static test equation (3) for this system is: 


{utl={Ft} 


where 

K33 = K^3p - g3 

F34 = ^ 34 p ■ ^ 3 , 94 

It can be seen that static test data is required on the first two bulkheads only 
if one of the first six rows are used to solve the first set of simultaneous 
equations. 


Test degree of freedom 



'"It 

^^2t 

"3t 

^"4t 

'"5t 

""6t 

'"Tt 

'"8t 

""9t 

'^lOt 

“lit 

“l2t 

1 

Kll 

Ki2 

0 

^14 

^15 

^16 

^17 

^18 

0 

Kl,10 

^1,11 

^1,12 

2 

^21 


^23 

^24 

^25 

^26 


K28 

0 

^2,10 

^2,11 

K 2 .I 2 

3 

0 

^32 

^33p 


^35 

^36 

0 

0 

K39 

^3,10 

^3,11 

^3,12 

4 

^41 

^42 

^43 

^44 

^45 


^47 


K49 

^4,10 

^4,11 

^4,12 

5 

^51 

^52 

^53 

^54 

^55 

^56 


K58 

K59 

^5,10 

^5,11 

^5,12 

6 

K61 

^62 

^63 

^64 

^65 


^67 

^68 

^69 

^6,10 

^6,11 

^6,12^ 


V 





/ 

V 





y 


First Bulkhead Second Bulkhead 

System Degree of Freedom System Degree of Freedom 

123456 123 45 6 
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Test data available . - The test data available on the payload consisted of a 
cantilever test, a bending test and a torque test. Any one of these tests is ade- 
quate provided sufficient data is measured in each test. The data needed and the 


data measured is a follows: 

Cantilever Bending Torque 

Test Test Test 



PS, ® P 


P & 6 measured 

Oj, O2 &62 also needed. 
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The major test displacements are shown on the previous page. However, all 
the degrees of freedom of the 1st and 2nd bulkhead having coupling terms in 
Figure 3 should be measured. The static test data available was not adequate 
because it only measured the vertical displacement of the first bulkhead (6i) 
in the cantilever test, the central deflection ( 6) in the beam bending test and 
the torsional rotation of the first bulkhead '(^i) in the torque test. More precise 
test data is required, but only on the first two bulkheads in this example. For 
the cantilever test and the bending test, row 3 would have to be the first solution 
to have the finite values (F or P) on the right hand side. If the torque static test 
IS used, row 4 would have to be used for the first solution. The degrees of 
freedom required for the measurements in each static test are: 


Cantilever test 
(Row 3) 

Bending test 
(Row 3) 

Torque test 
(Row 4) 


Bulkhead No. 1 
2, (3)^,4, 5, &6 

2, 3, 4, 5, & 6 

1, 2, 3, (4)^; 5, 

6 


Bulkhead No. 2 
3, 4, 5, & 6 

3, 4, 5, & 6 

1, 2, 3, 4, 5, & 6 
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c. Values that were measured in static test. 



The number of system measured eigenvalues and eigenvectors 
required depends on the row selected for the first solution. As previously 
stated the third row would be selected for the cantilever test or the bending 
test so as to have a force on the right hand side in the simultaneous equations 
solution. The fourth row would be used for the torque test. There are 
eleven unknown K terms in row 3 plus 9 unknown M terms or a total of 20 
unknowns. In row 4 there are 12 unknown M terms for a total of 24. 

If we used row 3 for the first set of equations to solve, the number of 
flexible modes required would be 13. There are 6 rigid body modes and 
one test equation. In examining the number of measured flexible modes 
available, it was found to be just short of the 13 that had all 6 degrees of 
freedom measured. The important point here is that for this 288 system 
degrees of freedom only 13 measured eignevalues and eigenvectors are 
required to solve the first set of simultaneous equations. Examining the 
coupled terms in the vehicle structure indicates that more measured 
eigenvalues may be required to solve the vehicle mass and stiffness matrix. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


1. A desynthesis procedure has been developed to combine static test 
data from the payload component with a few payload/ vehicle system eigen- 
values and mode shapes to derive the stiffness matrix and mass matrix 

of both components. The mass and stiffness matrix of each component is 
then used to extract the component eigenvalues and mode shapes. 

2. The desynthesis procedure was demonstrated on a two beam, 12 
degree of freedom system with nearly complete success. The last row 

of the payload component, and M^^ terms in equation (16) and the last 

row of the vehicle component K^2 12 ^12, 12 terms in equation (22) 

resulted in erroneous results as shown in Table 2. It is recommended 
that the order of row solution be revised so that the 5th row and 6th row 
solutions are interchanged and the 11th row and 12th row solutions are 
interchanged. The row solution order would then be 1, 2, 3, 4, 6, 5, 7, 

8, 9, 10, 12 and 11. Equations (15), (16), (21) and (22) would have to be revised 
and the computer code would have to be changed to reflect these alterations. 

3. It is recommended that the computer program be updated so that 
the variables can be interchanged on input to automate the following. 

. Change the order of row evaluation by an input variable. 

. By input variables, change the number of rigid body modes 
used, the number of flexible modes used, and the number 
of static test equations used in each set of row simultaneous 
equations. 

. Change the number of degrees of freedom of the system by 
input variables. 

. Change the location of coupled terms in each row by input 
variables 
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With these changes made the computer program would be automated to 
allow iterations On the solution of a given system and would have the 
capability to evaluate any system size and/or coupled degree of freedom 
location within the stiffness and mass matrix. 

4. It has been demonstrated that the accuracy of the solution is very 
sensitive to the accuracy of the static test data. It was demonstrated that 
in a simple beam bending test of the payload component, that a 1% error in 
the slope measurement or deflection in the center independently would 
cause considerable errors in M and K values. However, if they both had 

a 1% high measurement error, the correct M and K terms were derived. 
This desynthesis method has a large potential provided the accuracy that 
can be measured in the laboratory is sufficient. It is recommended that 
static test data be measured on a laboratory payload structure component 
and that this payload with a laboratory vehicle structure have the eigen- 
values and mode shape measured. This data would then be used in this 
desynthesis analysis to determine if the measurement accuracy of the test 
data is adequate. This laboratory structure could be a simple two beam 
model or the existing Langley Dynamic Research Model. 

5. The preliminary evaluation of the Z88 degree-of-f recdom Langley 
Dynamic Research Model consisting of a payload and a vehicle structure 
indicated that only 27 measured system eigenvalues and mode shapes and 
one static test of the payload should be sufficient to perform the desynthesis 
analysis. The payload static test could be a cantelever test, a simple 
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beam bending test or a cantelever torque test. Exact measurements 
of 6 degree of freedom displacements are required on two bulkheads at 
one end of the structure in each test. No measurements are required on 
the remaining 6 bulkheads. 

6. This desynthesis procedure should be applicable to a multicomponent 
system with a static test on one component. 


42 



TABLE 1 


EIGENVALUES AND EIGENVECTORS USED IN THE 
DESYNTHESIS ANALYSIS 

(VALUES FROM NASTRAN NORMALIZED TO ONE) 



i 

2 

3 

4 

1 

1,0 

1.0 

-. 9793407 

-. 8183043 

2 

0 

-1.0 

-. 0038746 

.6158865 

3 

1.0 

0 

-. 9271466 

. 0688562 

4 

0 

-1. 0 

. 1287581 

1. 0 

5 

1.0 

-1.0 

-. 7688752 

, 7703924 

6 

0 

-1. 0 

. 1225218 

. 2846892 

7 

1.0 

1. 0 

1. 0 

. 8224029 

8 

0 

-1, 0 

-. 0556343 

-. 7372945 

9 

1.0 

0 

. 9164300 

-. 0512733 

10 

0 

-1.0 

-. 1216034 

-. 9325306 

11 

1. 0 

-1.0 

. 7802890 

-. 8388912 

12 

__ 0 

-1.0 

-. 1181919 

-. 5815038 

EIGENVALUES 




0.0 


2 

X5 

1. 526401 (10) 

2 

*“2 

0. 0 



3. 078905 (10) 

2 

“3 

2. 155839 (10) 

5 

2 

U)7 

1. 367141 (10)' 

2 

0)4 

1. 674222 (10) 

6 

6 

00 

2. 625948 (10)' 


5 

6 

7 

8 

-. 4514553 

. 1480432 

-.2259740 

. 0368018 

1. 0 

-. 2254144 

1. 0 

-. 1144421 

. 2874545 

-. 1117015 

-. 0051818 

. 0058164 

-. 0258628 

. 0497821 

-. 8684909 

. 2216226 

-.4519291 

. 1266950 

. 2222345 

-.0226015 

-. 9194526 

. 0297036 

. 9141464 

. 0413581 

-. 1748622 

-.4540260 

-. 0435387 

-.2237451 

. 4297894 

1. 0 

, 2246833 

1.0 

. 1002866 

. 2895597 

. 0039955 

-. 0033242 

. 0165073 

-. 0248525 

-. 1292969 

-, 8428792 

. 1002866 

-. 4542416 

. 0396152 

.2217029 

-. 4588584 

-. 9219962 

. 2503231 

. 9470547 



TABLE 2 


ACCURACY OF STIFFNESS AND MASS TERMS 
GENERATED BY DESYNTHESIS PROCEDURE 


PAYLOAD COMPONENT 





KiZ 

Ki3 

^14 

^<i5 


Row 

1 

100. 0 

100. 0 

100. 0 

100. 0 


— 

Row 

2 


100. 0 

100. 0 

100. 0 

- 

- 

Row 

3 



100. 0 

- 

100.0 

100. 0 

Row 

4 




100. 0 

100. 0 

100. 0 

Row 

5 





100. 0 

100. 1 

Row 

6 

- 





99. 1 



M,i 

M.2 

M,3 

M^4 

Mi5 

M,6 

Row 

1 

'lOO. 0 

100. 0 

100. 0 

100. 0 

- 

_ 

Row 

2 


100. 0 

100. 0 

100. 0 

- 

- 

Row 

3 



100. 0 

- 

99. 9 

99. 9 

Row 

4 




100. 0 

99.8 

99.9 

Row 

5 





100. 0 

100. 1 

Row 

6 






84. 1 

J 


VEHICLE COMPONENT 




^^i7 

Ki8 

^19 

^ilO 

K.ll 

^12 

Row 

7 

100. 0 

100. 0 

100. 0 

100. 0 

_ 


Row 

8 


100. 0 

100. 0 

1 00. 0 

- 

- 

Row 

9 



100. 0 

- 

100. 0 

100. 0 

Row 

10 




100. 0 

100. 0 

100. 0 

Row 

11 





100. 0 

100. 1 

Row 

12 






100. 4 



- 





J 




m.h 

M,9 

^iio 

Mhi 

M ,, 
1 1 2 

Row 

7 

j 100. 0 

100. 0 

100. 0 

100. 0 

_ 

_ 

Row 

8 


1 00. 1 

99. 8 

99. 9 

- 

- 

Row 

9 



100. 1 

- 

99.9 

99. 9 

Row 

10 




100. 1 

100. 0 

100. 0 

Row 

11 





100. 0 

1 00. 0 

Row 

12 

L 





103.0 


INTERFACE SPRINGS 


= 100.0%: = 100.0%: ,, rr 99.4%: = 102.9% 

Nuniberh m matrix ri'f) i i>scnl tlio |K‘ri«'nt .iciui.iry compart'd to the values 
generated by the NASTRAN program 

(lyj, _ K or M (Desynthesis) ^ ^ qq 
K oi M (NASTRAN) 
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TABLE 3 


EFFECT OF STATIC TEST MEASUREMENT ERROR 
ON K AND M, SLOPE MEASURED 1% HIGH ON BOTH ENDS 

PAYLOAD COMPONENT 



'<ii 




^15 

^i6 


Row 1 

'l03. 1 

103. 1 

103. 1 

103. 1 

- 

- 


Row 2 


101. 4 

103. 0 

106. 7 

- 

- 


Row 3 



112.6 

- 

123. 3 

128.2 


Row 4 




107. 9 

119.4 

130. 1 


Row 5 





126. 9 

152. 8 


Row 6 





- 

265. y 



Mil 

Mi 2 

Mi3 

Mi 4 

M.5 

Mi6 


Row 1 

"l03. 1 

103. 1 

103. 1 

103. 1 




Row 2 


104. 3 

84.2 

94. 8 

- 

- 


Row 3 



97. 8 

- 

32. 2 

48. 7 


Row 4 




90. 3 

5. 3 

42. 0 


Row 5 





90.6 

144. 2 


Row 6 





- 

7676.- 


VEHICLE COMPONENT 






K,7 


^i9 

^10 

K 

ill 


Row 7 

103. 1 

103. 1 

103. 1 

103. 1 

- 

- 


Row 8 


101.6 

103. 1 

106. 3 

- 

- 


Row 9 



103.6 

- 

104.2 

106. 

0 

Row 10 




101.8 

102. 3 

104. 

2 

Row 11 





106. 6 

no. 

1 

Row 12 






226. 

8 


Mi7 

M,8 


^ilO 

Mill 

^il2 

Row 7 

103. 1 

103. 1 

103. 1 

103. 1 

- 

- 


Row 8 


111.4 

75. 7 

86.6 

- 

- 


Row 9 



106. 1 

- 

91. 3 

79. 6 

Row 10 




107. 8 

93. 2 

94. 

9 

Row 11 





76. 6 

30. 


Row 12 






13252. 


Numbers are the percentage related to the static test with no error 

K or M (Slope I7o High) 

% = X 100 

K or M (No Error) 


* Negative - 


indicates the values had opposite signs. 



TABLE 4 


EFFECT OF STATIC TEST MEASUREMENT ERROR 
ON K AND M, DEFLECTION MEASURED 1% HIGH AT CENTER 

PAYLOAD COMPONENT 



Kii 

K._ 

i2 

^i3 

Ki4 

Ki5 

Ki6 


Row 1 

" 96.2 

96.2 

96.2 

96.2 

- 

- 


Row 2 


97. 8 

96.2 

92.7 

- 

- 


Row 3 



88. 1 

- 

78.9 

74.3 


Row 4 




92. 3 

82.6 

73.4 


Row 5 





75. 1 

52. 0 


Row 6 






432.2 



Mil 


^.3 

i4 

M.(. 

i5 

^i6 


Row 1 

96.2 

96.2 

96.2 

91.2 

- 



Row 2 


95. 1 

114. 0 

103. 9 

- 

- 


Row 3 



101. 2 

- 

159. 8 

145. 5 


Row 4 




107.2 

180. 7 

149. 0 


Row 5 





108. 8 

63. 1 


Row 6 

- 





7191. 


VEHICLE COMPONENT 






Ki7 

^i8 


^ilO 


K.,_ 

il2 

Row 7 

96. 2 

96. 2 

96.2 

96.2 

- 

_ 


Row 8 


97.6 

96. 2 

93. 1 

- 

- 


Row 9 



95. 6 

- 

95. 1 

93. 

4 

Row 10 




97. 4 

96.9 

95. 

1 

Row 11 





92. 8 

89. 

4 

Row 12 






-17. 

¥ 


M_ 

i7 

^i8 

“,9 

“iio 


M 

il2 

Row 7 

^ 96.2 

96. 1 

96. 2 

96.2 

- 

- 


Row 8 


88. 3 

122. 0 

114. 5 

- 

- 


Row 9 



93. 3 

- 

107. 3 

118. 

3 

Row 10 




91.7 

105. 5 

103. 9 

Row 11 





122. 5 

169. 5 

Row 12 






-10035>! 


Numbers are the percentage related to the static test with no error 


% = 


K or M (Deflection 1% High) 
K or M (No Error) 


X 100 


* Negative 


indicates the values had opposite signs. 

4f> 



TABLE 5 


EFFECT OF STATIC TEST MEASUREMENT ERROR 
ON K AND M, END SLOPES AND CENTER DEFLECTION MEASURED 

1% HIGH 

PAYLOAD COMPONENT 



^ii 

KiZ 

Ki3 

Ki4 


Ki6 

Row 1 

■ 99.0 

99.0 

99.0 

99.0 

- 

- 

Row 2 
Row 3 
Row 4 
Row 5 
Row 6 


99. 0 

99.0 

99.0 

99.0 

99.0 

99.0 

99.0 

99.0 

99.0 
99.0 
99.0 
99. 0 


M.i 

Miz 


Mi4 

Mis 

Mi6 

Row 1 

r 99.0 

99.0 

99.0 

99.0 

- 


Row 2 
Row 3 
Row 4 
Row 5 
Row 6 

L. 

99. 0 

99.0 

99.0 

99.0 

99.0 

99.0 
99.0 
99. 0 

99.0 

99.0 

99.0 

99.0 


VEHICLE COMPONENT 




Ki7 

Ki8 

Ki9 

K,10 

Kill 

KilZ 

Row 

7 

^ 99.0 

99.0 

99.0 

99.0 



Row 

8 


99.0 

99.0 

99.0 

- 

- 

Row 

9 



99. 0 

- 

99. 0 

99. 0 

Row 

10 




99. 0 

99.0 

99.0 

Row 

11 





99. 0 

99.0 

Row 

IZ 






99.0 

•mJ 



M,7 


M-9 

M,i0 

^vll 



Row 7 

99.0 99.0 

99.0 

99. 0 

- 

- 

Row 8 

99. 0 

99. 0 

99. 0 

- 

- 

Row 9 


99. 0 

- 

99.0 

99.0 

Row 10 



99.0 

99.0 

99.0 

Row 11 




99.0 

99.0 

Row 12 





99.0 


Numbers are the percentage related to the static test with no error 


% 


K or M (Slope Sc Deflection 1% High) 
K or M (No Error) 


X 100 
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TABLE 6 


VEHICLE COMPONENT EIGENVALUES 
USING SCAMP 



Us 

ing De synthesis M and K 

Mode 

Using NASTRAN 

Analysis No. 1 

Analysis No. 2 

Analysis No. 4 

No. 

M and K 
(Hz) 

NO ERROR 
(Hz) ** 

Slope 
1% High 
(Hz) * 

Slope & Deflection 
1% High 
(Hz) 

3 

782. 5 

910.4 

- 

910.4 

4 

2449. 0 

2656.0 

- 

2656. 0 

5 

6124. 0 

6074.0 

- 

6074. 0 

6 

9783.0 

9254. 0 


9254. 0 


* Mass matrix not positive definite so SCAMP could not extract eigenvalues 

The discrepancy between NASTRAN and the desynthesis values is due to 
the fact that numerical problems were encountered in solving for Kj 2 12 
and Mj 2 12 (see Table 2). The Mj 2 j2 percent high, and the 

Ki 2 j 2 ^s 0.4 percent high. ^s discussed in Section IV, this dis- 

crepancy is probably due to the order of solution of K and M terms. 
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TABLE 7 


ACCURACY OF STIFFNESS AND MASS TERMS 
GENERATED BY DESYNTHESIS PROCEDURE 
(WITH ROW SOLUTION ORDER 1, 2, 3, 4, 6, 5, 7, 8, 9, 12, AND 11) 


PAYLOAD COMPONENT 




Kil 

K,2 

K.3 

K.,4 

K,5 

Ki6 

Row 

1 

r 

100. 0 

100. 0 

100. 0 

100.0 



Row 

2 


100. 0 

100.0 

100.0 

- 

- 

Row 

3 



100. 0 

- 

100. 0 

100.0 

Row 

4 




100.0 

100. 0 

100.0 

Row 

5 





100. 1 

100. 0 

Row 

6 






100. 1 



Mil 

M,2 

M,3 

M,4 

M,5 


Row 

1 

^100. 0 

100. 0 

100. 0 

100. 0 

_ 


Row 

2 


1 00. 0 

100.0 

100. 0 

- 

- 

Row 

3 



100. 0 

- 

99. 9 

99. 9 

Row 

4 




100.0 

99.8 

99.9 

Row 

5 





100. 1 

100. 0 

Row 

6 






100.^ 

VEHICLE COMPONENT 







Kj7 

K.8 

Ki9 

K,10 

^ ill 

K^12 



f 






Row 

7 

100. 0 

100. 0 

100. 0 

100. 0 

- 

Row 

8 


100. 0 

100. 0 

100. 0 

- 

- 

Row 

9 



100.0 

- 

100. 0 

100. 0 

Row 

10 




100. 0 

100. 0 

100. 0 

Row 

1 1 





100. 1 

100. 1 

Row 

12 

V. 





100. ,J 



M,7 



MilO 

Mill 

Mii2 

Row 

7 

^100. 0 

100. 0 

100.0 

100.0 

_ 

> 

Row 

8 


100. 1 

99. 8 

99. 9 

- 

- 

Row 

9 



100. 1 

- 

99. 9 

99.9 

Row 

10 




100. 1 

100. 0 

100. 0 

Row 

11 





100. 2 

100.2 

Row 

12i 






100. 7j 


INTERFACE SPRINGS 

Ki7 = 100.0%; K28 = 100.0%; ^ = 97.6%; 1%^ 12 = 100. 1%. 

Numbers in matrix represent the program accuracy compared to the values 
generated by the NAS'l’RAN program 

K or M (Desynthesis) 

= K or M (NASTPAN) 
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1 3 5 


4)2 4)4 4 ' 



l^ty^ rH 


1 = 3. 330 (10)-4 cm4 (8. 0 (10)-6 in4) 
E = 86185 N/cm2 (125000 psi) 
Density = 27. 68 g/cm3 (0. 1 lb/in3) 

A = . 06452 cm2 (. OlO in2) 

D = 2. 54 cm (1.0 in) 

El = 28. 7 N - cm2 (1. o lb - in 2) 


FIGURE 1.. TWO BEAM MODEL 
12 D. O.F. SYSTEM 
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Figure 1. (Continued) 


STIFFNESS MATRIX - 2 BEAM MODEL 
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Figure 1. (Concluded) 






MASS MATRIX - 2 
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FIGURE 2. - PAYLOAD - LANGLEY DYNAMIC RESEARCH MODEL 
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APPENDIX 


EFFECT OF SEQUENCE OF SIMULTANEOUS EQUATIONS 
ON ACCURACY OF STIFFNESSES AND MASSES 

As recommended in the second paragraph of the conclusions, the 
sequence of simultaneous equations used in the solution of equation (8) was 
revised so that stiffnesses and masses in the sixth matrix row were solved 
before those in the fifth row and, similarly, stiffnesses and masses in the 
twelfth row were solved before those in the eleventh row. Thus, the new 
row solution order used was 1, 2, 3, 4, 6, 5, 7, 8, 9, 10, 12, and 11. The 
original row order was 1 through 12 in numerical sequence- The results of 
the revised analysis are tabulated in Table 7 and should be compared with 
the original results in Table 2. By making this change, the inaccurate 
results for M^^, K|2^ 12^^^ ^12, 12 Table 2 were eliminated. All 

of the mass and stiffness terms of the payload component and the launch 
vehicle component were accurate to within one percent of the values from 
the NASTRAN solution. The interface spring stiffness term had the 

largest error, namely 97.6 percent of the NASTRAN value. 
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